Generations of orthogonal surface coordinates

Two generation methods were developed for three dimensional flows where the computational domain normal to the surface is small. With this restriction the coordinate system requires orthogonality only at the body surface. The first method uses the orthogonal condition in finite-difference form to determine the surface coordinates with the metric coefficients and curvature of the coordinate lines calculated numerically. The second method obtains analytical expressions for the metric coefficients and for the curvature of the coordinate lines

Re-entrant ferroelectricity in liquid crystals

The ferroelectric (Sm C$^*$) -- antiferroelectric (Sm C$^*_A$) -- reentrant ferroelectric (re Sm C$^*$) phase temperature sequence was observed for system with competing synclinic - anticlinic interactions. The basic properties of this system are as follows (1) the Sm C$^*$ phase is metastable in temperature range of the Sm C$^*_A$ stability (2) the double inversions of the helix handedness at Sm C$^*$ -- Sm C$^*_A$ and Sm C$^*_A$% -- re-Sm C$^*$ phase transitions were found (3) the threshold electric field that is necessary to induce synclinic ordering in the Sm C$^*_A$ phase decreases near both Sm C$^*_A$ -- Sm C$^*$ and Sm C$^*_A$ -- re-Sm C$^*$ phase boundaries, and it has maximum in the middle of the Sm C$^*_A$ stability region. All these properties are properly described by simple Landau model that accounts for nearest neighboring layer steric interactions and quadrupolar ordering only.Comment: 10 pages, 5 figures, submitted to PR

Generations of orthogonal surface coordinates

An orthogonal surface-oriented coordinate system has been developed for three-dimensional flows where the computational domain normal to the surface is small. With this restriction the coordinate system requires orthogonality only at the body surface. This coordinate system has been utilized in boundary layer flows and for the hypersonic viscous shock-layer problem. Two methods have been developed for generating the surface coordinates. 9 figures. (RWR

Chemically Reacting Viscous Flow Program for Multi-Component Gas Mixtures.

A general computer program was developed for solving the laminar boundary layer equations with a finite-difference method. The governing equations are solved in an uncoupled manner in order that a gas mixture with a large number of chemical species can be readily handled. The program has been written with various options to provide a flexibility that allows a variety of problems to be solved with only a change in the input data

Correlated defects, metal-insulator transition, and magnetic order in ferromagnetic semiconductors

The effect of disorder on transport and magnetization in ferromagnetic III-V semiconductors, in particular (Ga,Mn)As, is studied theoretically. We show that Coulomb-induced correlations of the defect positions are crucial for the transport and magnetic properties of these highly compensated materials. We employ Monte Carlo simulations to obtain the correlated defect distributions. Exact diagonalization gives reasonable results for the spectrum of valence-band holes and the metal-insulator transition only for correlated disorder. Finally, we show that the mean-field magnetization also depends crucially on defect correlations.Comment: 4 pages RevTeX4, 5 figures include

Coupled Flow Field Simulations of Charring Ablators with Nonequilibrium Surface Chemistry

This paper describes the coupling of a Navier-Stokes solver to a material response code to simulate nonequilibrium gas-surface interactions. The Navier-Stokes solver used in this study is LeMANS, which is a three-dimensional computational fluid dynamics code that can simulate hypersonic reacting flows including thermo-chemical nonequilibrium effects. The material response code employed in this study is MOPAR, which uses the one-dimensional control volume nite-element method to model heat conduction and pyrolysis gas behavior. This coupling is demonstrated using a test case based on the Stardust sample return capsule. Coupled simulations are performed at three different trajectory conditions. The effects of the pyrolysis gas chemistry are evaluated by assuming that the gas is either in chemical equilibrium or composed entirely of non-reacting phenol. The results show that the non-reacting pyrolysis gas assumption produces higher convective heat fluxes, surface temperatures, and mass blowing rates. These effects are mainly due to the composition of the pyrolysis gas. The additional species produced by the pyrolysis gas in the chemical equilibrium case react with oxygen and nitrogen atoms in the gas-phase. This results in fewer atoms participating in the exothermic surface reactions, which reduces the heat transfer to the vehicle

Continuum breakdown in hypersonic viscous flows

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76654/1/AIAA-2002-651-530.pd

Computational Modeling of Gas-Surface Interactions for High-Enthalpy Reacting Flows

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106443/1/AIAA2013-187.pd

On Validating an Astrophysical Simulation Code

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.Comment: 76 pages, 26 figures (3 color), Accepted for publication in the ApJ

Conceptual Analysis of Electron Transpiration Cooling for the Leading Edges of Hypersonic Vehicles

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/140441/1/6.2014-2674.pd